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Abstract 



We study B — * rj'X s within the framework of the Standard Model. Sev- 
eral mechanisms such as b — > rj'sg through the QCD anomaly, and b — ► rfs 
and S — > r/sg arising from four-quark operators are treated simultaneously. 
Using QCD equations of motion, we relate the effective Hamiltonian for the 
first mechanism to that for the latter two. By incorporating next-to-leading- 
logarithmic(NLL) contributions, the first mechanism is shown to give a signif- 
icant branching ratio for B — * rj'X s , while the other two mechanisms account 
for about 15% of the experimental value. The Standard Model prediction for 
B — * i]'X s is consistent with the CLEO data. 
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The recent observation of B — > rj'K QXj| and B — > r)'X s decays with high momentum 7/ 
mesons has stimulated many theoretical activities |p|-|lO|]. One of the mechanisms proposed 
to account for this decay is £» — ^ sg* — > sgrj' [^ffl where the rj' meson is produced via 
the anomalous rf — g — g coupling. According to a previous analysis J4j, this mechanism 
within the Standard Model(SM) can only account for 1/3 of the measured branching ratio: 



B{B -> rfX a 



6.2 ± 1.6(stat) ± 1.3(syst)+^(bkg 



x 10" 



with 2.0 < p v > < 2.7 GeV. 



There are also other calculations of B — > rj'X s based on four-quark operators of the effective 
weak-Hamiltonian [§,0. These contributions to the branching ratio, typically 10~ 4 , are 
also too small to account for B — > i]'X s , although the four-quark-operator contribution is 
capable of explaining the branching ratio for the exclusive B — ► rfK decays These 
results have inspired proposals for an enhanced b —> sg and other mechanisms arising from 
physics beyond the Standard Model in order to see if new physics should play any 

role in B —* r]'X s , one has to have a better understanding on the SM prediction. In this 
letter, we carry out a careful analysis on B — > f]'X s in the SM using next-to-leading effective 
Hamiltonian and consider several mechanisms simultaneously. 

We have observed that all earlier calculations on b — > sgrj' were either based upon one- 
loop result [|J which neglects the running of QCD renormalization -scale from My/ to M b or 
only taking into account part of the running effect ||. Since the short-distance QCD effect 
is generally significant in weak decays, it is therefore crucial to compute b — > sgrj' using 
the effective Hamiltonian approach. As will be shown later, the process b — > sgrj' alone 
contribute significantly to B — > rj'X s while contributions from b —>■ rj's and B — > rj'sq are 
suppressed. 

The effective Hamiltonian 



H] for the B — > rj'X s decay is given by: 



H eff (AB = 1)= % £ V fb V* s {CMO f M + CMO f M) 

-KVHCECiiriOiQi) + c 8 (/i)o 8 (/i))], 

i=3 



(1) 



with |2 



0\ — {-Sij 'j)v -a(J : jbi)v -A, J 2 — {Sifi)v-A{fjbj)v-A 

3 = {sibi) V -Aj2(QjQj)v-A, 4 = {sib^v-AY^iW^v-A 
q 

5 = (Sibi)v-A ^.(QnQn)v+A, Oh = (Sib 

Os 



f 



-•<>>< I \ r -A 22(<ljQj)v+A, 
q 



iUj)V-A ^2{qjQi)v+A 

q 



^L- Si a» v (m s P L + m b P R )T«b 3 G%, 



47T 



(2) 



where V ± A = 1 ±75. In the above, we have dropped O7 since its contribution is negligible. 
For numerical analyses, we use the scheme-independent Wilson coefficients discussed in Ref. 
130 ■ For m t = 175 GeV, a s (m|) = 0.118 and jj, = m b = 5 GeV, we have pj) 



Ci 



-0.313, C 2 = 1.150, C, = 0.017, C 4 



-0.037, C s = 0.010, C fi 



-0.045, 



(3) 



At the NLL level, the effective Hamiltonian is modified by one-loop matrix elements which 
effectively change C$(//)(i = 3, • • ■ , 6) into Cj(//) + Ci(q 2 , fi) with 
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C 4 (q 2 ,fi) = C- 6 (g 2 ,/i) = -3C 3 (q 2 ,^ = -3C 5 (q 2 ,n) = -P s (q 2 ,fi), (4) 



where 



with 



P s (g 2 ,/i) = ^-CM (y + G(m 2 ,g 2 ,/i) ) . (5) 



G(m 2 , q 2 , p) = 4 / x(l - x) log ( ^ ^ x)q * ) dx. (6) 



The coefficient C§ is equal to —0.144 at /i = 5 GeV [|TTJ, and m c is taken to be 1.4 GeV. 

Before we discuss the dominant b — > sgrj' process, let us first work out the four-quark- 
operator contribution to B — > i]'X s using the above effective Hamiltonian. We follow 



the approach of Ref. which uses factorization approximation to estimate various 



hadronic matrix elements. The four-quark operators can induce three types of processes 
represented by 1) < ri^qT^B >< X s \sT[q\0 >, 2) < T]'\qT 2 q\0 >< X s \sTb\B >, and 3) 
< r]'X s \sT 3 q\0 >< 0\qT' 3 b\B >. Here denotes appropriate gamma matirces. The con- 
tribution from 1) gives a "three-body" type of decay, B — > rj'sq. The contribution from 2) 
gives a "two-body" type of decay b — > sr]'. The contribution from 3) is the annihilation 
type which is relatively suppressed and will be neglected. Note that there are inteferences 
between 1) and 2), so they must be coherently added together ||. 

Several decay constants and form factors needed in the calculations are listed below: 

< 0|u 7/ ,7 5 m|V >=< 0\dj^ 5 d\r]' >= iftp* 

2 

f 1 / "^Tii / 

< 0|s7 M 7 5 s|r/ >= if*,pl , < 0|s7 5 s|7? >= i(f% - f^')^, 
f% = cos 9 1 + -j=fa sin 9 8 ), = -^=(h cos 9 1 - V2f 8 sin 8 ) , 

— 777 /~?^ 777^ 

< V '\ui»b\B- >=< r/|J 7 ^° >= Fftof+tf) + (F* - F^) ^— 5 

= sin^ 8 + cosM^). (7) 

Fot the rf — rj mixing associated with decay constants above, we have used the two-angle 



?2 2 

V 
-Vim, 



-parametrization. The numerical values of various parameters are obtained from Ref. |T6 
with fx = 157 MeV, / 8 = 168 MeV, and the mixing angles 9 X = -9.1°, 9 8 = -22.1°. 
For the mixing angle associated with form factors, we use the one-angle parametrization 
with 9 = —15.4° ||16]| , since these form factors were calculated in that formulation ||[15 



In the latter discussion of b — > sgrj', we shall use the same parametrization in order to 
compare our results with those of earlier works For form factors, we assume that 

pBrn _ pBrjs _ pB-K clip i e anc [ monopole q 2 dependence for F\ and F , respectively. 
We used the running mass m s ~ 120 MeV at fi = 2.5 GeV and F B?r = 0.33 following Ref. 

The branching ratios of the above processes also depend on two less well-determined 
KM matrix elements, V ts and V u b. The dependences on V ts arise from the penguin- diagram 
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contributions while the dependences on V u b and its phase 7 occur through the tree-diagram 



contributions. We will use 7 = 64° obtained from Ref. [17], \V ts \ ~ Wcb\ = 0.038 and 
|Kt6|/|K*| = 0-08 for an illustration. We find that, for [l = 5 GeV, the branching ratio in 
the signal region p v > > 2.0 GeV {m x < 2.35 GeV) is 

B{b -> r)'X s ) ss 1.0 x 10" 4 . (8) 

The branching ratio can reach 2 x 10 -4 if all parameters take values in favour of B — > T]'X S . 
Clearly the mechanism by four-quark operator is not sufficient to explain the observed 
B — > r]'X s branching ratio. 

We now turn to the major mechanism for B — > rj'X s : b — ► rj'sg through the QCD 
anomaly. To see how the effective Hamiltonian in Eq. ([1]) can be applied to calculate this 
process, we rearrange part of the effective Hamiltonian such that 

QOi = (C 3 + ^)0 3 + (C 5 + ^)0 5 - 2(C 4 - C 6 )O a + 2(C 4 + C 6 )O v , (9) 



where 

A = s Tm (1 - 7 5 )T a 6^g 7 M 75T a g, O v = ^(1 - l5 )T a bJ2^T a q. (10) 

Since the light-quark bilinear in Oy carries the quantum number of a gluon, one expects |3j 
Oy give contribution to the b — > sg* form factors. In fact, by applying the QCD equation 
of motion : D U G£" = g s T,qi^T a q, we have O v = (l/g^sj^l - ^)T a bD u G^ 0. In this 
form, Oy is easily seen to give rise to b —>■ sg* vertex. Let us write the effective b — ► sg* 
vertex as 

If = -^VZV a ^(AFiStf7* - i ^) Ln ~ ^F 2 m b sa^RT a b). (11) 

In the above, we define the form factors AF\ and F2 according to the convention in Ref. [f§] . 
Inferring from Eq. (|9]), we arrive at 

AF 1 = —(CM + CM), F 2 = -2CM (12) 



We note that our relative sign between AFi and F 2 agree with those in Ref. [§,0], and shall 
result in a destructive interference for the rate of b — ► sgrj'. We stress that this relative 
sign is fixed by treating the sign of Os and the convention of QCD covariant derivative 
consistently fl9| . To ensure the sign, we also check against the result by Simma and Wyler 
20| on b — » sg* form factors. An agreement on sign is found. Finally, we remark that, 



at the NLL level, AF% should be corrected by one-loop matrix elements. The dominant 
contribution arises from the operator 2 where its charm-quark-pair meets to form a gluon. 
In fact, this contribution, denoted as AF\ for convenience, has been shown in Eqs. (4)-(6), 
namely AF t = £(C* 4 (g 2 , //) + C 6 (q 2 , //)). 

To proceed further, we recall the distribution of the b(p) — > s(p') +g(k) + rj'(k') branching 
ratio H: 
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FIG. 1. The distribution of B(b — > s + g + rf) as a function of the recoil mass mx- 
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lAFxl'co + i?e(AFiF 2 *)— + \AF 2 \ 2 ^ 

y y 2 



(13) 



where a g (p) = V^Vf«s(/^) l^frf is t ne strength of 77' — g — g vertex: a g cos Oe^a^k 13 with q 
and k the momenta of two gluons; x = (p' + k) 2 /ml and y = (k + k') 2 /m 2 ; Co, c\ and C2 are 
functions of x and y as given by: 



c 

Cl 
C-2 



-2x 2 y +(l-y)(y- x')(2x + y- x')] /2, 
(l-?/)(l/-x') 2 , 

2x 2 t/ 2 - (1 - y){y - x'){2xy -y + x')} /2, 



(14) 



with x' = m 2 //m 2 ; and the — r] mixing angle is taken to be —15.4° as noted earlier. 
Finally, in obtaining the normalization factor: 0.2, we have taken into account the one-loop 
QCD correction [21] to the semi-leptonic b —>■ c decay for consistency. 

In previous one-loop calculations without QCD corrections, it was found Ai*\ w —5 and 
F 2 w 0.2 [|,g. In our approach, we obtain AF 1 = -4.86 and F 2 = 0.288 from Eqs. (|) and 
(|I"2|). However, AFi is enhanced significantly by the matrix-element correction AFi(q 2 ,fi). 
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The latter quantity develops an imaginary part as q 2 passes the charm-pair threshold, and 
the magnitude of its real part also becomes maximal at this threshold. From Eqs. (3), (4) 
and (5), one finds Re(AF 1 (4m 2 , fi)) = —2.58 at \i = 5 GeV. Including the contribution by 
AFx(g 2 ,/i) with \i = 5 GeV, and using Eq. (|i~3"D, we find B(b — > sgrj') = 5.6 x 10~ 4 with the 



cut mx = y(k + p') 2 < 2.35 GeV imposed in the CLEO measurement 0. This branching 
ratio is consistent with CLEO's measurement on the B —>■ rj'X s branching ratio [TJ]. Without 
the kinematic cut, we obtain B(b — > sgrj') = 1.0 x 10~ 3 , which is much larger than 4.3 x 10~ 4 
calculated previously ||. We also obtain the spectrum dBib — > sgr]')/dm x as depicted in 
Fig. 1. The peak of the spectrum corresponds to mx ~ 2.4 GeV. 

It is interesting to note that the CLEO analysis [§] indicates that, without the anomaly- 
induced contribution, the recoil-mass(mx) spectrum of B — > i]'X s can not be well reproduced 
even if the four-quark operator contributions are normalized to fit the branching ratio of the 
process. On the other hand, if b — » sg* — > sgrj' dominates the contributions to B — > i]'X s , 
as shown here, the mx spectrum can be fitted better as shown in Fig. 2 of Ref. [[|. It 
is also interesting to remark that although the four-quark operator contributions can not 
fit the branching ratio nor the spectrum, it does play a role in producing a small peak 
in the spectrum, which corresponds to the B — > rj'K mode. Specifically, the B — > rj'K 
mode is accounted for by the b — > srj' type of decays discussed earlier. Based on results 
obtained so far, one concludes that the Standard Model is not in conflict the experimental 
data on B — > f]'X s . It can produce not only the branching ratio for B — > r]'X s but also the 
recoil-mass spectrum when contributions from the anomaly mechanism and the four-quark 
operators are properly treated. 

Up to this point, a s (/x) of the rf — g — g vertex has been treated as a constant independent 
of invariant-masses of the gluons, and \i is set to be 5 GeV. In practice, a g {p) should behave 
like a form-factor which becomes suppressed as the gluons attached to it go farther off- 
shell However, it remains unclear how much the form-factor suppression might be. 
It is possible that the branching ratio we just obtained gets reduced significantly by the 
form-factor effect in rf — g — g vertex. Should a large form-factor suppression occur, the 
additional contribution from b — > rj's and B — ► rj'sq discussed earlier would become crucial. 
We however like to stress that our estimate of b — > sgrj' with a s evaluated at \i = 5 GeV 
is conservative. To illustrate this, let us compare branching ratios for b — > sgrj' obtained at 
fi = 5 GeV and /x = 2.5 GeV respectively. In NDR scheme |22], branching ratios at the above 
two scales with the cut mx < 2.35 GeV are 4.9 x 10~ 4 and 9.1 x 10 -4 respectively. One can 
clearly see the significant scale-dependence! With the enhancement resulting from lowering 
the renormalization scale, there seems to be some room for the form-factor suppression in 
the attempt of explaining B —>■ 1]'X S by b —>■ sgrj' |23| . 



It should be noted that the above scale-dependence is solely due to the coupling constant 
cx s (n) appearing in the rf — g — g vertex. In fact, the b —>■ sg* vertex is rather insensitive 
to the renormalization scale. Indeed, from Eq. (|ll|), we compute in the NDR scheme the 
scale-dependence of g s ■ (AFj + AFi(q 2 )). We find that, as \x decreases from 5 GeV to 2.5 
GeV, the peak value of the above quantity increases by only 10%. Therefore, to stablize the 
scale-dependence, one should include corrections beyond those which simply renormalize the 
b — > sg* vertex. We shall leave this to a future investigation. 

It is instructive to compare our results with those of Refs. |3]|§]. With the kinematic 
cut, our numerical result for Bib — > sgrj') is only slightly smaller than the branching ratio, 
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8.2 x 10~ 4 , reported in Ref. 0, where the a s (fi) coupling of rf — g — g vertex is evaluated at 
H ~ 1 GeV, and AF\ receives only short-distance contributions from the Wilson coefficients 
C4 and Cq. Although we have a much smaller a s , which is evaluatd at fi = 5 GeV, and 
the interference of AFi and F 2 is destructive Q rather than constructive [|J, there exists a 
compensating enhancement in AF± due to one-loop matrix elements. The branching ratio 
in Ref. || is 2 — 3 times smaller than ours since it is given by a AF± smaller than ours 
but comparable to that of Ref. ||. Concerning the relative importance of AFi and F2, we 
find that Ai*\ alone gives B(b — > sgrj') = 6.5 x with the kinematic cut mx < 2.35 
GeV. Hence the inclusion of F 2 lowers down the branching ratio by only 14%. Such a small 
interference effect is quite distinct from results of Refs. []3],|| where 20% — 50% of interference 
effects are found. We attribute this to the enhancement of Aiq in our calculation. 

Before closing we would like to comment on the branching ratio for B — > r]X s . It is 
interesting to note that the width of b — > rjsg is suppressed by tan 2 9 compared to that of 
b — > rj'sg. Taking 9 = —15.4°, we obtain BiB — > f]X s ) k4x 10 -5 . The contribution from 
the four-quark operator can be larger. Depending on the choice of parameters, we find that 
B(B — > r]X s ) is in the range of (6 ~ 10) x 10~ 5 . 

In conclusion, we have calculated the branching ratio of b — > sgrj' by including the NLL 
correction to the b — > sg* vertex. By assuming a low-energy 7/ — g — g vertex, and cutting 
the recoil-mass mx at 2.35 GeV, we obtained B{b — > sgrj') = (5 — 9) x 10~ 4 depending on 
the choice of the QCD renormalization-scale. Although the form-factor suppression in the 
1]'— g—g vertex is anticipated, it remains possible that the anomaly-induced process b sgr]' 
could account for the CLEO measurement on BiB — > r]'X s ). For the four-quark operator 
contribution, we obtain BiB — > r]'X s ) k 1 x 10~ 4 . This accounts for roughly 15% of the 
experimental central-value and can reach 30% if favourable parameters are used. Finally, 
combining contributions from the anomaly-mechanism and the four-quark operators, the 
entire range of B — > r]'X s spectrum can be well reproduced. 
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